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ABSTRACT 



The delivery truck problem is one In which a truck is 
loaded with m packages, one package to be delivered to each 
of m destinations. The amount of fuel consumed by the truck 
is directly dependent upon the current total v/eight of the 
truck, which includes both the v/eight of the packages and 
the amount of fuel remaining in the tank. The problem is 
to determine a sequence in which to deliver all m packages 
which will minimize total fuel consumption. A branch and 
bound algorithm for obtaining optimal solutions to the 
delivery truck problem is presented, along with several 
sample problems with their solutions. A brief report 
of computational experience is included. 
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I. INTRODUCTION 



The delivery truck problem is one in which a truck is 
loaded with m packages, one package to be delivered to each 
of m destinations. The amount of fuel consumed by the 
truck is directly dependent upon the current total vjeight 
of the truck, which includes both the v/eight of the packages 
and the amount of fuel remaining in the tank. The problem 
is to determine a sequence in which to deliver all m package 
which v;ill minimize total fuel consumption. 

This problem is closely related to the traveling 
salesman problem. It differs in that the route is not 
closed and that fuel consumed on each leg of the journey 
is not only dependent upon the length of the leg, but also 
the current vjeight of the truck. 

This thesis is the presentation of an algorithm for 
obtaining optimal solutions to the delivery truck problem. 
The methodology of branch and bound is used to implicitly 
enumerate all possible solutions. Several example problems 
with solutions are presented followed by a brief report 
of computational experience for problems with m ranging 
from five to twenty. 



II . RELATED RESEARCH 



The delivery truck problem is clearly a problem in 
the area of discrete mathematical optimization. Reference 
8 is a current summary of results applicable to such 
problems . 

Current research on the ordinary traveling salesman 
problem is summarized in Ref. 2 , in which, after an analysis 
of computational results, it is recommended that dynamic 
programming [Ref. 5] be used for problems with thirteen 
cities or less, and that branch and bound [Ref. 9] be used 
for the rest. 

An algorithm for obtaining optimal solutions to the 
modified traveling salesman problem in which no return to 
the starting city is required is presented in Ref. ?• 

In an extension of the traveling salesman problem called 
"The Delivery Problem" [Ref. H2 more than one salesman is 
available to visit the predetermined set of cities. 

Another extension of the traveling salesman problem, 
sometimes referred to as the "Milk Bottle Problem", [Ref. 
is that in which the milkman, (salesman), is to deliver 
one or more bottles of milk to a known set of houses, 
(cities). He must walk and wishes to minimize the amount 
of v;ork expended carrying the milk bottles. The problem 
is very closely related to the delivery truck problem, 
the difference being that no account is taken for fuel 
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used along the way. A branch and bound algorithm for the 
milk bottle problem is presented in Ref. 3. This algorithm 
is a direct extension of the work of Little, et al., [Ref. 
6 ]. 

The delivery truck problem can be formulated as an 
investigation problem. Current results in Investigation 
Theory are contained in Ref. 1. The algorithm presented 
in this thesis uses an approach similar to that of the 
branch and bound algorithm presented in Ref. 1. 
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III. MATHEMATICAL DESCRIPTION 



The following notation vjlll be used. Let 



(d 



IJ 



) 



7T 



f 

7r 



f 

7T 



i 






= a matrix of distances betv^een destination 
1 and destination j, i = l,2,...,n, 

= (iT^, it 2 > • . • , iT^) be a permutation of 
the Integers l,2,...,n representing an 
order in which destinations l,2,...,n 
are visited. 

= the amount of fuel required if path tt 
is followed. 

= the amount of fuel required if destinations 
1 T^, tt 2 > are visited in that 

ohder. Clearly f = f 

7T TT^ 

= fuel required to travel from destination 
IT. to destination it. given that 1 destina- 
tions have already been visited and f 
fuel has been consumed. 



Let f =0. Then 



IT. 



O 



j=i-l 

f 2 F (j>F )• 

""l j = 0 ""j+l 



The delivery truck problem is to select that path tt which 
minimizes f . 

TT 
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The algorithm contained herein assumes it is known 
which destination v;lll be the first visited, and the problem 
starts with the truck initially located there. If the 
first destination is not known, the truck can be assumed 
to be located at a dummy location added to the given set 
of destinations, with travel allowed to any valid destination 
at zero fuel consumption. 

Given that the truck is initially located at destination 
one, there are (n-1)! possible solutions. 

In different variations of the delivery truck problem, 
different fuel consumption functions F are appropriate. In 
this thesis, F is taken to have the follovjlng form: 



F 



IT 



j 



(i 






f 




c^(c2 





)d 



IT . IT 
1 
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where 



c^.= an arbitrarily chosen constant 

C 2 = the fully loaded weight of the truck 

c^ = the weight of package k, k = l,...,n. 



Note that if f =0 for all i, the problem reduces 
to the milk bottle problem of Ref. 3* If> In addition, 
c^ = 0, i = l,...,n, the problem reduces to the modified 
traveling salesman problem in which no return to the 
starting city is required. 
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F 




Fi^:;ure 3. 1 

For the experimental testing discussed in Section IV(E) 
the constants chosen were c^ = .004, = 8 OOO and 

c^ = 270 , for 1 = i,2,...,n. The resulting function is 
used to Illustrate variation in fuel consumption due to 
truck loading as shown in Figure 3.1* Line a represents 
consumption with 1=0 and f = 0, line b for 1 = 5, f = 0, 
and line c for i = 5 and f = 500. 
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IV. THE ALGORITHM 



This section begins by describing how the branch and 
bound method is applied to the delivery truck problem, 
follo\ved by a general statement of the algorithm. (A flow 
diagram of the algorithm is contained in Appendix A and 
the computer program immediately aftervrards . ) Several 
sample problems are used to illustrate the output of the 
algorithm. A brief summary of computational experience 
is presented. 

A. APPLICATION OF BRANCH AND BOUND TO THE DELIVERY 
TRUCK PROBLEM 

The branch and bound method consists of selectively 
partitioning the set of all feasible solutions and computing 
bounds on the objective function for each element of the 
partition. This process is continued until one element 
of the partition containing a single solution is obtained 
for which the associated bound is at least as good as that 
of any other element. 

A branching tree is initiated and extended such that 
each node at the end of a branch of the tree corresponds 
to a set of solutions. The set of all such nodes specifies 
a partition of the solution space. The bounds associated 
with each element of the partition, and therefore with 
each node at the end of a branch of the tree, represent a 
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lower bound on the amount of fuel used if any solution 

t 

contained in that element of the partition is followed. 

Associated with every node of the tree is a destination. 
The initial node, representing all solutions, corresponds 
to destination number one, the Initial location of the 
truck. Branching corresponds to selecting a node from 
which to branch and in specifying a destination for use 
in extending the path specified by that branch. The branch 
corresponding to node one naturally is selected initially 
for extension. This branch is extended to the right and 
the left by adding two nev; nodes to the tree. See Figure 
4.1. The node labeled k represents all solutions in 




Figure 4.1 Initial branching 

which the truck starts at destination one and proceeds first 
to destination k. The node labeled -k represents all solutions 
in which the truck does not go to destination k first. 
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Associated v/lth each node at the end of a branch of 



the tree is a path, starting with destination one, and 
including each destination number corresponding to a 
positively labeled node in the order in v;hlch they appear 
in the branch. 

When a branch of the tree is selected for extension, 
the path specified by that branch is called the current 
path. 

The bound on each ending node with a positive label 
applies to all solutions starting vjlth the path specified 
by that node. The bound on those with negative labels 
corresponds to all solutions starting with the path specified 
v/hich do not include as their next element, all negative 
labels of nodes at the end of that branch. In illustration 
of this, consider Figure ^.2. The bound on node labeled 5 




Figure 4.2 Sample tree 
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applies to all solutions starting v;ith elements 1, 3 , 5 , in 
that order. The bound on the node labeled -5 applies to 
all solutions starting with elements 1, 3 > and not having 
7 or 5 as the next element. 

After a branch is selected for extension, it is necessary 
to determine vihlch destinations are eligible for use in 
extending the current path. This is done by determining 
all destinations not in the current path and not prohibited 
due to negatively labeled nodes at the end of the branch. 

The order of an element of the partition reduces to 
one v;hen there are no destinations eligible for use in 
extending the branch corresponding to that element. V/hen 
such an element is located with a bound less than or equal 
to the bound on all others, the path corresponding to the 
sole solution contained in that element is optimal. 

B. GENERAL STATEMENT OF THE ALGORITHM 

Specific methods of branching and bounding are ignored 
temporarily while the general structure of the algorithm 
is presented. 

Algorithm: 

1. Create the initial node corresponding to destination 
number one. Compute a bound for this node. Go to 2. 

2. Select a branch to extend. This branch determines the 
current path. Go to 3. 

3. Determine the set E of destinations eligible for use 

in extending the current path. Put into E all destinations 
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not in the current path. Remove from E all indices corres- 
ponding to negatively labeled nodes at the end of the branch 
being extended. If |E| = 0, stop. If not, go to 4. 

4. Select from E a destination k for use in extending the 
current path. Extend the branch to the left and right 
creating two new nodes labeled k and -k. Go to 5* 

5. Compute bounds for branches ending with nodes labeled 
k and -k. Go to 2. 

Upon termination the current path is optimal. 

C. BRANCHING AND BOUNDING RULES 
1 . Branching 

Branching is a two part operation, consisting first 
of selecting a branch to extend, and next of selecting a 
destination for use in extending that branch. It has 
proved to be advantageous in branch and bound to use as 
the first of these steps the selection of that branch 
whose associated bound is smallest and specifying tie 
breaking rules in the event they occur. This rule is used 
in the algorithm. 

The second part of the branching operation selects 
a destination for use in extending the current path. The 
algorithm selects that eligible destination which is closest 
to the last element in the current path, or, in other words, 
the unvlsited destination closest to the truck's current 



position . 



2 . Bounding 



The bound associated with each branch specifies a 
lovjer bound on the amount of fuel used If the path specified 
by that branch Is follovjed. In Illustration of the method 
used, suppose branching Is being carried out from a node 
labeled j to destination k, and that there are r destinations 
In the current path. The exact amount of fuel required to 
follow the current path Is f . , and the fuel required to 

. J 

travel from destination j to destination k Is F, (r, f.). 

k ’ j 

The lovjer bound on the node labeled k Is 

^ fj) f L 

where L Is a lower bound on the fuel required to visit the 
remaining n-r-1 destinations. 

L Is obtained by considering the upper triangle of 
the reduced matrix ) where 1 and j are not destinations 

In the current path, selecting the smallest element In 
row k and the remaining n-r-2 smallest elements, ordering 
these distances In non-decreasing order, and computing the 
amount of fuel required to traverse these distances In 
the prescribed sequence. 

Only the upper triangle need be considered since 
(d . ) Is symmetric and travel between destination 1 and 

1 J 

destination J Is carried out at most once. 

It can be seen that L Is a valid lower bound by 
noting that after leaving destination k, the truck must 
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traverse an additional n-r-1 legs in its required journey. 
Clearly the shortest distance it will be required to haul 
n-r-1 packages is the smallest element in rov; k of the 
reduced matrix. A lov:er bound on the next leg is clearly 
the fuel required to haul n-r-2 packages the smallest 
distance in the reduced matrix, excluding the element used 
for the first leg. This procedure is continued until all 
n-1 legs of the journey are accounted for. 

When the bound is computed for the node labeled -k, 
the procedure is the same except that selection of a distance 
for the first additional leg after leaving destination j 
is not restricted to come from rov/ k of the reduced matrix. 

D. PROOF OP OPTIMALITY 

Let S be the set of all solutions to the problem, i.e., 

S = {it|tt is a permutation of integers l,2,...,n} , 
and let P be a partition of S such that 




Let b. be a lov;er bound on all solutions contained in 
1 

P^, i.e.. 



b. < f for all TT in P 

1 — IT 



16 



In Step 2 of the algorithm a branch of the tree is 
selected for extension. Selecting a branch for extension 
corresponds to selecting an element of the current partition 
for refinement. If element P. is selected and IP.| = 1, 

0 ‘o'’ 

no further refinement is possible, and in Step 3, |e| = 0. 

At this point all n destinations are contained in the 

current path, and bound b. is no longer a lower bound, but 

0 

exact. Note that the method of selection used for P., as 
described in Section IV(C)1, is to select partition element 
j such that 

b. = min b, for all P, in P. 
j k k 

Hence, P^ has an exact bound which is less than or equal 
to the lov;er bounds associated with all elements of P. 

The single solution it in P^ is thus a solution which 
minimizes f . 

TI 

E. SAMPLE PROBLEMS AND SOLUTIONS 

The following are sample problems v;ith optimal solutions. 
The X and y coordinates of the destinations were selected 
from a uniform distribution over the interval (0, 10). 
Coordinates were adjusted slightly when necessary to 
insure that each pair of destinations v;as at least a unit 
distance apart. The truck is assumed to be initially 
located at destination one. The path shov/n is the one which 
visits each of the remaining destinations and uses the 
minimum amount of fuel. 
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Sannle Problorr 1 
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So.rv.^le Problen 2 
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Sa’'r''?le 



Problem 3 
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F. COMPUTATIONAL EXPERIENCE 



Computational experience vjith the algorithm is limited, 
but an investigation was made into the v;ay the time required 
to solve an n destination problem varies with n. Seventeen 
sample problems were developed as described in the preceeding 
section; five with n equal to five, six w'ith n equal to 
ten, and six w’ith n equal to fifteen. The average times 
required to solve each set of problems, along v;ith the 
average number of bounds computed, is shown in Table ^.1. 



n 


average 
time ( sec . ) 


average § 
of bounds 


5 


.38 


5.1 


10 


1.6 


145.0 


15 


99.7 


2525.0 



Table ^.1 



Although these results are sparse, this data suggests 
that time required increases exponentially with n, as seen 
in Figure 4.9. The circles show the data points of Table 
4.1. The dots show times required to solve sample problems 
having n equal to 20. 
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Figure 4.9 



22 



APPENDIX A 



FLOW DIAGRAM OP COMPUTER PROGRAM 
The following notation, not previously defined, is used 



LN 


= the last node created 


HN ( I ) 


= destination number associated v;ith node I, 
for 1=1, 2 , . . . , LN 


PRED(I) 


= the node immediately node I in the branching 
tree, for 1=2, . . . , LN 


BFN 


= node selected for extension 


BTH 


= destination selected for use in extending 
the current path 


BB(I) 


= bound on node I 



B(I) = BB(I) is node I is at the end of a branch in the 





tree, and equals M otherwise, where M is 
large . 


NPL 


= number of packages currently delivered 


MN 


= destination at which the truck is currently 
located 


P(X,Y,Z) 


= fuel required to travel distance Z if X 
packages have currently been delivered and 
Y fuel consumed doing so. 
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BRANCH ANC 60LND ALGORITHM FOR THE DELIVERY TRUCK PROBLEM 



C IMENS ICN X(l5),Y(15),NE(15),NEE(15),B(4CCCJtBB(A0C0), 
INFM 400CJ ,NPRED ( ^»OOC) ,NZ(^OOC) , ZZ (4000) , PCS ( L5 ) , 

IN ROW ( 15) ,NCOL{ 15 ) , DN 15 , 15 ) , 00 ( 1 5 , 15 ) 

C 

READ N, THE NUMBER OF DESTINATIONS 
C 

RE AD (5, SCO N 
90C FORMAT (12) 

C 

READ (X,Y), THE CGCRDINATES FOR EACH CESTINATICN 
C 

DC 11 I=1,N 
READ(5,9C1 ) X( I) ,Y(I) 

9C 1 FCRMAK 2F5.2) 

NRITE(6,333) I,X(I),Y(I) 

333 FCRMATC I = ',I4,' X=',F5.2,' Y=*,F5.2) 

11 CONTINUE 

CO 750 1=1, N 
DO 750 J = 1 ,N 
IF( I .Eg .J.) GO TO 751 

C ( I , J) = ( (X ( I )-X( J J )**2+( YI I )-Y( J ) . 5 

DD( I , J) =D( I , J) 

GO TO 75C 
751 C ( I , J)=99999. 

DC( I , J) =99999. 

75C CONTINUE 
LN-1 
B( 1 )=C. 

EE ( 1)=0 .0 
MN = 1 

NFN( 1)=1 
NF= 1 

NFPEDI 1 )=1 
DC 10 I=1,N 
N E E ( I ) = 1 
1C NE(I)=1 
NBFN=1 
NETh=l 
NFN = 2 
NFL=1 
NNN=N-1 
l^T = 2200 . 

101 IF(I\FN(NBFN) .LT.O) GO TO 120 
DM IN=99999. 

DNEXT = 999.99. 

DC 701 J=2,N 

IF(NE( J) .£0.0) GO TO 701 
IF(D(NP , J) .GE .DNEXT) GO TO 701 
IF(C(NP, J) .GE.DMIN) GO TO 7C2 
DNEXT=DMI N 
NEXT=NBTH 
CM IN=D( NP, J ) 

NBTH=J 
GO TO 7C1 
7C2 CN£XT=D (NP,J ) 

NEXT=J 

7C1 CONTINUE 
LN=LN+1 
INTP=LN/1CC 
NR = LN-INTP'UOO 
IF(NR.NE.O) GO TO 760 
759 IF (BB(NFN) .LT ,999.0, ) GO TO 760 
NFN=NFN+1 
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GC TO 759 
76C CCMiKUE 

N1-N'(LN} =KBTH 
KFREC(LN)=NBFN 

E ( LKj =b (MM + .004*(8000 .-NPL*270 .-e(MN) ) «CMIK 
END=d( LN) 

KK=N-NPL-i 
K ~ C 

7C2 SPALL=999S9. 

CC 704 1=2, NNN 
L = I + 1 

DC 704 J=L,N 

IT- {NEE( Ih.EQ.O) GO TO 704 
IF(NEE( JJ.EG.O) GC TO 704 
IF( SMALL.LE.D( I , J1 ) GC TO 7C4 
SFALL=U(I, JJ 
I I=I 
JJ = J 

704 CCKTINUE 
K = K + 1 

PCS(K)=SPJ\LL 
NPCW(K}=I I 
NCCL( K) =JJ 
D( I I , J J } = 99999. 

IF (K.LT .KK ) GO TO 702 
DC 705 1=1 ,KK 
I 1 I=NROR( I ) 

J JJ=NCGL( I ) 

C ( I I 1, J JJI =CD(I I I , JJJ ) 

7C5 CONTINUE 
K=1 

JX=NPL+1 

ice FF=.G04'^( 8CCC.-JX’I'27C.-BNDJ*P0S (K) 

ENC=BN04FF 
K = K+1 
JX = JX+ I 

I F (K.L E .KK ) GO TO 7C6 
BE ( LN) =END 

707 LN=LN+1 
INTP=LN/ 100 

NR = LN-I NTP=!=100 
IF(NR.NE.O) GO TO 77C 

769 IF(Bb(l\FNJ .LT .9990. ) GO TO 770 
NFN=NFN+1 

GC TO 7E9 

770 CCKTINUE 
KFK(LN) =-NBTH 
KPREC{LNJ=KBFK 

EKC=B( MK )+ .004* (80 00. -NPL -270.-6 (MN) )*DNEXT 
K — C 

709 SKALL=99999. 

CC 708 1=2, NNN 
L=I + 1 

CC 708 J=L.N 

(f (hEEi D:. EC.O) GO TO 708 
IF (NEE( J) . EQ.O) GC TO 708 
IF(SMALL.LE.D( I ,J) ) GO TO 708 
SRALL=D(I,J) 

I I = I 
J J=J 

708 CCKTINUE 
K = K + 1 

PCS(K)=SMALL 
KRCV«(K) = M 
KCCL(K) =JJ 
C( I I ,JJ )=99999. 

IF(K.LT .KK) GO TO 709 
DC 710 1=1 ,KK 
I I I=NROK(.I ) 

JJJ=NCOL( I ) 

C(III,JJJ)=CD(I II,JJJ) 

71C CCKTINUE 
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K=1 

JX=KFL+1 

711 FF=.004=!^( 8C00.- JX*270. -BND ) *POS ( K) 
EI\C = BNO + FF 
K = K + 1 
JX=JX+ 1 



IF(K.LE.KK) GO TO 711 
B0{ LN) -6KD 
NZ(LN1=KEXT 
ZZ(LM=CKEXT 
BB (^BFNJ =99999. 

GC TO 120 
120 NBTF!=NZ(i\BFN1 
LF = 1,K+1 
IKTP=LN/LCC 
KP=LN-I KTP^lOO 
IF (KR.NE.jO) go TC 780 

779 IF(83(i\:FN) .LT.9990. J GO TO 780 
KFK=NFN+1 

GC TO 779 

780 CONTINUE 
(\FK(LN1 =NBTFi 
NPRED(LN).-NBFN 

E (LN)=b (MM + .004=5= ( 8000. -NPL*270.-E(MN) ) =!=ZZ( NBFN) 
ENC = B( LM 
KK=N-NPL-1 
K = 0 

712 SFALL=999,99. 

DC 713 1=2, NNN 
L=I + 1 

DC 713 J=L,N 

IF(NEE( I l-EQ.O) GO TO 713 
1F(NEE( J) .EG.OJ GO TO 712 
1F(SMALI.LE.C(I , Jl) GC TO 713 
SMALL=0( I, J) 



I I = I 
JJ = J 

712 CONTINUE 
K = K + 1 

FCS(K)=SMALL 
NROW(K) = M 
NCCL(K)=JJ 
C ( I I ,JJ )=99999. 

IF(K.LT.KKJ GO TO 712 
CC 714 1=1, KK 
I I I=NROW( I ) 

J JJ=NCCL( I } 

C(1I I, J JJ1=DD(I I I, JJJ ) 

71^ CONTINUE 
K=1 

JX=NFL+1 

715 FF=.004*(8C0C.-JX*270.-BND)*P0S (K) 
ENC=BND+FF 

K = K + 1 
JX=?JX+1 

IF(K.LE.KK) GO TO 715 
EE(LN)=BND 
NE(NBTH 1=0 
CNEXT=99999 . 

DC 716 J=2,N 

IF(NE( J) .EC.O) GC TC 716 
1F(C(NP, J) .GE.ONEXT) GO TO 716 
DNEXT = D (NP , J) 

NEXT=J 

716 CONTINUE 
GC TO 7C7 

120 XMINB=9999. 

CC 131 I = i\F(\,LN 

IF(BB( I ) .GE.XMINB) bO TC 131 

XM IN6=BB( I 1 

N6FN=I 

131 CONTINUE 
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DC 140 1=2, K 
K' E E ( n = 1 

14 0 NE ( n = i 

l^^ = ^3F h 

141 fviP=KHN(KN) 

IF (KP.GT .0 ) GO TO 150 

KC=-KP 

KE(NO}=C 

PK=fvPREC (MM 

GC TO 141 

15C RPL=1 
PP = MN 

151 IF (MM. EG. 1) GC TO 160 
K=KHN(MMI 

IF(K.LT .0) GC TO 152 

MEE(K) =C 

NE(K)=0 

MFL=NPL+1 

152 MM=KPPEC(MM} 

GC TO 151 

160 IFdvPL.EG.M GO TO 170 
^ELF=0 

DC 161 1=2, N 
NELH=NELHi-KE ( 1 ) 

1F(NELF.GT.U GO TO 101 
IF(NE( II .EG.O) GC TC 161 
LAST1=I 

161 CCNTINUE 
LK=LK+1 
IMP=LN/1CC 
RP=LN-I MTP^lOO 
IF(KR.NE.O) GO TC 730 

72S IF(BB(NFN) ,LT.9990) GO TO 730 
N FK = NFN-»1 
GC TO 729 

730 CCKTINUE 

M-M( LNJ =LAST1 
K'PRED( LM =KBFM 

e(LIM) = B (MN ) + .004«( 8000 .-NPL*270 .-B ( MN I I =i=D ( KP , LAST 1 1 

eB(LN)=B(LM 

BB(NBFN)=99999. 

GC TO 130 

170 RFIT£(6,949) 

949 FCRMAT(‘1*,‘ OPTIMAL SEQUENCE, IK REVERSE ORDER',///) 
KX=KBFK 

800 KK = KHK'(KBFM 

UFITE(6,950) NN 

950 FORMAT! 110,/) 

IF(NBFN.EQ.l) GO TO 830 

831 KBFK=NPRtD (KBFN) 

IF(NhN( KBFM .GT. 0) GO TO 800 
GC TO 831 

830 CCKTINUE 

KRITE(6,951) LN 

951 FORMAT!///, • NUMBER OF edUNDS C OMPLTED= ‘ , I 6 ) 
RRITE(6,952) B(NX) 

952 FCRMAT(F12.2) 
kRITE(6,960) NFN 

960 FORMAT! ' FN=' ,110) 

STCP 

EKC 
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